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Differential distributions for heavy quark production depend on the heavy quark mass and other 
momentum scales, which can yield additional large logarithms and inhibit accurate predictions. 
Logarithms involving the heavy quark mass can be summed in heavy quark parton distribution 
functions in the ACOT factorization scheme. A second class of logarithms involving the heavy- 
quark transverse momentum can be summed using an extension of Collins-Soper-Sterman (CSS) 
formalism. We perform a systematic summation of logarithms of both types, thereby obtaining an 
accurate description of heavy-quark differential distributions at all energies. Our method essentially 
combines the ACOT and CSS approaches. As an example, we present angular distributions for 
bottom quarks produced in parity-conserving events at large momentum transfers Q at the ep 
collider HERA. 
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INTRODUCTION 



In recent years, significant attention was dedicated to 
exploring properties of heavy-flavor hadrons produced in 
lepton-nucleon deep inelastic scattering (DIS). On the 
experimental side, the Hadron-Electron Ring Accelerator 
(HERA) at DESY has generated a large amount of data 
on the production of charmed |jl|, ||, [| ||, |^] and bottom 
mesons M, ffl, |l pi O. At present energies (of order 300 
GeV in the ep center-of-mass frame) , a substantial charm 
production cross section is observed in a wide range of 
Bjorken x and photon virtualities Q 2 , and charm quarks 
contribute up to 30% to the DIS structure functions. 

On the theory side, Perturbative Quantum Chromo- 
dynamics (PQCD) provides a natural framework for the 
description of heavy-flavor production. Due to the large 
masses M of the charm and bottom quarks (M 2 
Aq CD ), the renormalization scale can be always cho- 
sen in a region where the effective QCD coupling as is 
small. Despite the smallness of as, perturbative calcu- 
lations in the presence of heavy flavors are not without 
intricacies. In particular, care in the choice of a factor- 
ization scheme is essential for the efficient separation of 
the short- and long-distance contributions to the heavy- 
quark cross section. This choice depends on the value 
of Q as compared to the heavy quark mass M. The key 
issue here is whether, for a given renormalization and fac- 
torization scale hf ~ Q, the heavy quarks of the N-th 
flavor are treated as partons in the incoming proton, i.e., 
whether one calculates the QCD beta-function using N 
active quark flavors and introduces a parton distribution 
function (PDF) for the A^-th flavor. A related, but sepa- 
rate, issue is whether the mass of the heavy quark can be 
neglected in the hard cross section without ruining the 
accuracy of the calculation. 



Currently, several factorization schemes are available 
that provide different approaches to the treatment of 
these issues. Among the mass-retaining factorization 
schemes, we would like to single out the fixed flavor 
number factorization scheme (FFN scheme), which in- 
cludes the heavy-quark contributions exclusively in the 
hard cross section |l2] , |l3] , |l4[ [T5| , [l6| |; and massive 
variable flavor number schemes (VFN schemes), which 
introduce the PDFs for the heavy quarks and change the 
number of active flavors by one unit when a heavy quark 



threshold is crossed ll 18, 19, GO 21, 22, E3 



Further details on these schemes can be found later in 
the paper. Here we would like to point out that, were 
the calculation done to all orders of as, the FFN and 
massive VFN schemes would be exactly equivalent. How- 
ever, in a finite-order calculation the perturbative series 
in one scheme may converge faster than that in the other 
scheme. In particular, the FFN scheme presents the most 
economic way to organize the perturbative calculation 
near the heavy quark threshold, i.e., when Q 2 m M 2 . 
At the same time, it becomes inappropriate at Q 2 ^> M 2 
due to powers of large logarithms In (Q 2 /M 2 ) in the hard 
cross section. In the VFN schemes, these logarithms 
are summed through all orders in the heavy-quark PDF 
with the help of the Dokshitzer-Gribov-Lipatov-Altarelli- 
Parisi (DGLAP) equation [§e| |?[ §§|; hence the pertur- 
bative convergence in the high-energy limit is preserved. 
In their turn, the VFN schemes may converge slower at 
Q 2 ?s M 2 , mostly because of the violation of energy con- 
servation in the heavy-quark PDF's in that region. Re- 
cently an optimal VFN scheme was proposed that com- 
pensates for this effect po| . 

In this paper, we would like to concentrate on the anal- 
ysis of semi-inclusive differential distributions (i.e., dis- 
tributions depending on additional kinematical variables 
besides x and Q). We will argue that finite-order calcu- 



lations in any factorization scheme do not satisfactorily 
describe such distributions due to additional large loga- 
rithms besides the logarithms ln(Q 2 /M 2 ). To obtain sta- 
ble predictions, all-order summation of these extra loga- 
rithmic terms is necessary. 

The extra logarithms are of the form 
(a n s /q%) ln m (g|/Q 2 ), < m < 2n - 1, where q T = Pt/z 
denotes the transverse momentum pr of the heavy 
hadron in the 7*p center-of-mass (cm.) reference frame 
rescaled by the variable z = (pa ■ Ph)/(pa • q)- Here 
p A , q* 1 , and p^ H are the momenta of the initial-state 
proton, virtual photon, and heavy hadron, respectively. 
Our definitions for the j*p cm. frame and hadron 
momenta are illustrated by Fig. [I]. The resummation of 
these logarithms is needed when the final-state hadron 
escapes in the current fragmentation region (i.e., close 
the direction of the virtual photon in the j*p cm. 
frame, where the rate is the largest). In the current 
fragmentation region, the ratio q\jQ 2 is small; therefore, 
the terms \n m (q^/Q 2 ) compensate for the smallness 
of as at each order of the perturbative expansion. If 
hadronic masses are neglected, such logarithms can be 
summed through all orders in the impact parameter 
space resummation formalism [3l[ [$|[ |3^| , which 
was originally introduced to describe angular correla- 
tions in e + e~ hadroproduction [Q |35|] and transverse 
momentum distributions in the Drell-Yan process p^. 1 
Here the impact parameter b is conjugate to qr- The 
results of Refs. | p0[ [3l| , [32|, ^| are immediately valid 
for semi-inclusive DIS (SIDIS) production of light 
hadrons (it, K, ...) at Q of a few GeV or higher, and for 
semi-inclusive heavy quark production at Q 2 M . 
To describe heavy-flavor production at Q 2 ~ M 2 , the 
massless QT-resummation formalism must be extended 
to include the dependence on the heavy-quark mass M. 

In this paper, we perform such extension in the Aivazis- 
Collins-Olness-Tung (ACOT) massive VFN scheme JTp] 
with the optimized treatment of the threshold region P9[ . 
We adopt a "bottom-up" approach to the development of 
such mass-dependent resummation. 2 We start by sepa- 
rately reviewing the massive VFN scheme in the inclu- 
sive DIS and qr resummation in the massless SIDIS. We 
then discuss a combination of these two frameworks in 
a joint resummation of the logarithms ln(Q 2 /M 2 ) and 
ln(g^/(5 2 ). As a result, we obtain a unified descrip- 
tion of fully differential heavy-hadron distributions at all 
Q 2 above the heavy quark threshold. It is well known 
that the finite-order calculation does not satisfactorily 



The similarity between the multiple partem radiation in semi- 
inclusive DIS and the other two processes was known for a long 
time; see, for instance, an early paper |37| . 
2 An alternative "top-down" approach will require the analysis of 
leading regions in the high-energy limit and derivation of the evo- 
lution equations that retain terms with positive powers of M/Q. 
Such analysis could involve methods similar to those discussed 
in Ref. 
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Figure 1: The parity-conserving semi- inclusive production 
e+p — > e + H + X of heavy hadrons in the 7*p cm. reference 
frame. The resummation effects considered here are impor- 
tant in the current fragmentation region 6h — > 0, i.e., when 
the final-state heavy quark h closely follows the direction of 
its escape in the O(ag) flavor-excitation process 7* + h — > h. 



treat the current fragmentation region for any choice 
of the factorization scheme. In contrast, the proposed 
massive extension of the gT-resummation accurately de- 
scribes the current fragmentation region in the whole 
range Q 2 > M 2 . 

The present study is interesting for two phenomeno- 
logical reasons. Firstly, the quality of the differential 
data will improve greatly within the next few years. By 
2006, the upgraded collider HERA will accumulate an in- 
tegrated luminosity of 1 fb _1 |59|, i.e., more than eight 
times the final integrated luminosity from its previous 
runs. Studies of the heavy quarks in DIS are also en- 
visioned at the proposed high-luminosity Electron Ion 
Collider @ and THERA Q. Eventually these exper- 
iments will present detailed distributions both at small 
(Q 2 « M 2 ) and large (Q 2 ^> M 2 ) momentum transfers. 

Secondly, the knowledge of the differential distribu- 
tions is essential for the accurate reconstruction of in- 
clusive observables, such as the charm component of the 
structure function F 2 (x,Q 2 ). At HERA, 40-60% of the 
charm production events occur outside the detector ac- 
ceptance region, notably at small transverse momenta of 
the heavy quarks. To determine F£(x, Q 2 ), those events 
should be reconstructed with the help of some theoretical 
model, which so far was the 0(a 2 s ) calculation in the FFN 
scheme ||l5| , fl2"| , |43| , pi[ incorporated in a massless par- 
ton showering generator. As mentioned above, for inclu- 
sive observables the FFN scheme works the best not far 
from the threshold and becomes unstable at Q 2 3> M 2 , 
where the VFN scheme is more appropriate. In more de- 
tail, the VFN scheme describes F£(x, Q 2 ) better than the 
FFN scheme when Q 2 exceeds 20 (GeV/c) 2 , i.e., roughly 
when Q 2 /M 2 > 10 [|o). The transition to the VFN 
scheme occurs faster at smaller x, where the cm. energy 
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of the 7*p collision is much larger than M. For bottom 
quark production, the estimate Q 2 /M 2 > 10 corresponds 
to Q 2 > 200 (GeV/c) 2 . The VFN calculation can also be 
extended down to the mass threshold to uniformly de- 
scribe the whole range Q 2 > M 2 . Since the proposed re- 
summation calculation is formulated in the VFN scheme, 
it provides a better alternative to the finite-order calcu- 
lation in the FFN scheme due to its correct treatment of 
differential distributions at all values of Q 2 . 

As an example, we apply the developed method to the 
leading-order flavor-creation and flavor-excitation pro- 
cesses in the production of bottom mesons at HERA. We 
find that the resummed cross section for this process can 
be described purely by means of perturbation theory due 
to the large mass of the bottom quark. Our predictions 
can be tested in the next few years once the integrated 
luminosity at HERA approaches 1 fb _1 . Essentially the 
same method can be applied to charm production. In 
that case, however, the resummed cross section is sensi- 
tive to the nonperturbative large-6 contributions due to 
the smaller mass of the charm quarks, and the analysis 
is more involved. Since the goal of this paper is to dis- 
cuss the basic principles of the massive g/r-resummation, 
we leave the study of the charm production and other 
phenomenological aspects for future publications. 

The paper has the following structure. Section || re- 
views the application of the ACOT factorization scheme 
[|l8| and its simplified version jlT], to the description of 
the inclusive DIS structure functions. Section |l| recaps 
basic features of the 6-space resummation formalism in 
massless SIDIS. Section discusses modifications in the 
resummed cross section to incorporate the dependence on 
the heavy-quark mass M. In Section 0, we present a de- 
tailed calculation of the mass-dependent resummed cross 
sections in the leading-order flavor-creation and flavor- 
excitation processes. Section [v| presents numerical re- 
sults for polar angle distributions in the production of 
bottom quarks at HERA. Appendix [A| contains details on 
the calculation of the 0(as) mass-dependent part of the 
resummed cross section. In Appendix we present ex- 
plicit expressions for the 0(as) finite-order contributions 
from the photon-gluon channel. Finally, Appendix |c] dis- 
cusses in detail the optimization of the ACOT scheme 
when it is applied to the differential distributions in the 
vicinity of the threshold region. 



II. OVERVIEW OF THE FACTORIZATION 
SCHEME 

A. Factorization in the presence of heavy quarks 

In this Section, we discuss the application of 
the Aivazis-Collins-Olness-Tung (ACOT) factorization 
scheme Jl8) to inclusive DIS observables, for which this 
scheme yields accurate predictions both at asymptoti- 
cally high energies and near the heavy-quark threshold. 

In the inclusive DIS, the factorization in the presence 



of heavy flavors is established by a factorization theo- 
rem jLTjl, which we review under a simplifying assumption 
that only one heavy flavor h with the mass M is present. 
Let A denote the incident hadron. According to the theo- 
rem, the contribution F/,m(^, Q 2 ) of h to a DIS structure 
function F(x, Q 2 ) (where F(x, Q 2 ) is one of the functions 
Fi(x,Q 2 ) or Fi(x, Q 2 ) in parity-conserving DIS) can be 
written as 




Here the summation over the internal index a includes 
both light partons (gluons G and light quarks), as well 
as the heavy quark h. This representation is accurate up 
to the non-factorizable terms that do not depend on M 
and can be ignored when Q 3> Aqc d ■ The non- vanishing 
term on the r.h.s. is written as a convolution integral 
of parton distribution functions f a /A (£j {M-f/ to <}}) an d 
coefficient functions C h / a (x a /ti, Vf/Q, M/Q). The con- 
volution is realized over the hadron light-cone momen- 
tum fraction £ carried by the parton a. The coefficient 
function depends on the flavor-dependent "scaling vari- 
able" Xa discussed below. The parton distributions and 
coefficient functions are separated by an arbitrary fac- 
torization scale (Xf such that f a /A depends only on fip 
and quark masses {m q } = m u , mj, ...,M; and Ch/ a de- 
pends only on fip, M, and Q. As a result of this sepa- 
ration, all logarithmic terms Og ln fe ([ip/m q ) with light- 
quark masses are included in the PDF's, where they are 
summed through all orders using the DGLAP equation. 
Note that in the massless approximation such logarithms 
appear in the guise of collinear poles 1/e in the pro- 
cedure of dimensional regularization. The logarithms 
ln fe (/j,f/M) with the heavy- quark mass M are included 
either in Ch/ a or depending on the factorization 

scheme in use. 

In the reference frame where the momentum of the 
incident hadron A in the light-cone coordinates is 

(where p ± = (p° ±p 3 ) /v2), the quark PDF can be de- 
fined in terms of the quark field operators ip q (x) as |^5| 

/«4{£=})-rr/£***- 

x v y ' J spin color 

x(pA|^ 9 (0,y~,0 T ) 

xPcxp j-ig J dz~j2f + (0, z~,0 T )| 

x^ q (0)\ PA ). (2) 

Here Pexpj...} is the path-ordered exponential of the 
gluon field srf v {x) in the gauge rj ■ srf = 0. The r.h.s. 
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is averaged over the spin and color of A and summed 
over the spin and color of q. A similar definition exists 
for the gluon PDF. The dependence of / q /a(£> iVf/"^}) 
on fip is induced in the process of renormalization of 
ultraviolet (UV) singularities that appear in the bilo- 
cal operator on the r.h.s. of Eq. (||). In general, the 
PDF is a nonperturbative object; however, it can be 
calculated in PQCD when fiF 3> Aqcd, and the inci- 
dent hadron A is replaced by a parton. This feature 
opens the door for the calculation of F h / A (x, Q 2 ) for any 
hadron A through the conventional sequence of calcu- 
lating C/>/ a (x Q /£, hf/Q, M/Q) in parton-level DIS and 
convolving it with the phenomenological parameteriza- 
tion of the nonperturbative PDF f a /A(£,, {HF/m q }). In 
the inclusive DIS, it is convenient to choose \if ~ Q to 
avoid the appearance of the large logarithm In^p/Q) in 

C h/a (Xa/Z,HF/Q,M/Q). 

The factorized representation (|l|) is valid in all fac- 
torization schemes. The specific factorization scheme is 
determined by (a) the procedure for the renormalization 
of the UV singularities; and (b) the prescription for keep- 
ing or discarding terms with positive powers of M/ Q in 
the coefficient function Ch/ a . The choice (a) determines 

if the logarithms hn. k (^,F /M) are resummed in the heavy- 
flavor PDF or not. With respect to each of two issues, the 
choice can be done independently. For instance, the MS 
factorization scheme uses the dimensional regularization 
to handle the UV singularities, but does not uniquely de- 
termine the choice (b). Hence, it is not necessary in this 
scheme to always neglect M in the coefficient function 
and expose the heavy-quark mass singularities as poles 
in the dimensional regularization. 

The ACOT scheme belongs to the class of the vari- 
able flavor number (VFN) factorization schemes |46| that 
change the renormalization prescription when \ip crosses 
a threshold value fXthr- It is convenient to choose fithr 
for the flavor h to be equal to M, since the logarithms 
ln fc (fip/M) vanish at that point. If [If < M, all graphs 
with internal heavy-quark lines are renormalized by zero- 
momentum subtraction. If [ip > M, these graphs are 
renormalized in the MS scheme. The masses of the light 
quarks are neglected everywhere, and graphs with only 
light parton lines are always renormalized in the MS 
scheme. 

The physical picture behind the ACOT prescription is 
simple: the heavy quark is excluded as a constituent of 
the hadron for sufficiently low energy (an TV — 1 flavor 
subscheme), but the heavy quark is included as a con- 
stituent for sufficiently high energies (an TV flavor sub- 
scheme). The renormalization by zero- momentum sub- 
traction below the threshold leads to the explicit decou- 
pling of the heavy-quark contributions from light parton 
lines. As one consequence of the decoupling, all pertur- 
bative components of the heavy-quark PDF vanish at 
fiF < M, so that a nonzero heavy-quark PDF may ap- 
pear only through nonperturbative channels, such as the 
"intrinsic heavy quark mechanism" |Q. Since the size 
of such nonperturbative contributions remains uncertain, 



they are not considered in this study. At [If > M, a 
non-zero heavy-quark PDF fhlA is introduced, which is 
evolved together with the rest of the PDFs with the help 
of the mass- independent MS splitting kernels. The ini- 
tial condition for fh/A{£., A*f) is obtained by matching the 
factorization subschemes at fiF = M. At order as, this 
condition is trivial: 



fh/A(t,l*F = M) = 0. 



(3) 



At higher orders, the initial value of fh/A^i Mf) is given 
by a superposition of light-flavor PDF's ^0|. A simple 
illustration of these issues is given in Appendix [a|. 

The ACOT scheme possesses another important prop- 
erty: the coefficient function Ch/ a in this scheme has a 
finite limit as Q — > oo, which coincides with the expres- 
sion for the coefficient function obtained in the massless 
MS scheme with TV active flavors. This happens be- 
cause the mass-dependent terms in C^u contain only 
positive powers of M/Q, while the quasi-collinear log- 
arithms 1ii(/zf/M) are resummed in A*f). As a 
consequence of the introduction of fh/A(£i A*f), the coeffi- 
cient function C h / a includes subprocesses of three classes: 

• flavor excitation, where the parton a is a heavy 
quark; 

• gluon flavor creation, where a is a gluon; 

• and light-quark flavor creation, where a is a light 
quark. 

In contrast, in the FFN scheme @ |]J [l| [l| 

only the flavor-creation processes are present. The 
lowest-order diagrams for each class are shown in Fig. ||. 

The subsequent parts of the paper consider the pro- 
cesses shown in Figs, ^a and ||b. Note that we count the 
order of diagrams according to the explicit power of as 
in the coefficient function, i.e., 0{a s ) in Fig. ||a, 0{a s ) 
in Fig. ||b, and 0(a s ) in Fig. ||c. This counting does 
not apply to the whole structure function F h / A (x, Q 2 ) in 
Eq. (Q) when the heavy-quark PDF is itself suppressed by 
as/7r near the mass threshold p8| , |is| , pof . In that region, 
an O(og) flavor- excitation contribution has roughly the 
same order of magnitude as the 0(a 7 g + ) flavor- creation 
contribution. We return to this issue in the discussion 
of numerical results in Section (v|, where we interpret 
the combination of the 0(a s ) flavor- excitation contri- 
bution (Fig. ||a) and 0(a s ) flavor-creation contribution 
(Fig. ||b) as a first approximation at Q f» M. 



B. Simplified ACOT Formalism 

Of several available versions of the ACOT scheme, 
our calculation utilizes its modification advocated by 
Collins @, which we identify as the Simplified ACOT 
(S-ACOT) formalism |2^]. It has the advantage of being 
easy to state and of allowing relatively simple calcula- 
tions. This simplicity could be crucial for implementing 
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(a) 



(b) 



(c) 



Figure 2: Basic subprocesses in the ACOT scheme: (a) flavor 
excitation 7* + h — > h at 0(a%); (b) gluon flavor creation 
(photon-gluon fusion) 7* + G — > /i + ft, at C(a s ); (c) light- 
quark flavor creation 7* + q — + (7* + G) + q — > (ft + ft) + <? 
at C(q|). The thick and thin solid lines correspond to the 
heavy quark ft and light quarks q = u,d, s, respectively. 



the massive VFN prescription at the next-to-leading or- 
der in the global analysis of parton distributions. 
In brief, this prescription is stated as follows. 

Simplified ACOT (S-ACOT) prescription: 
Set M to zero in the calculation of the coef- 
ficient functions C^/a for the incoming heavy 



quarks: that is, 



s~t I A,n t-"r 
^h/h I — , -^T, — 



Q 



M 
Q 



a 



h/h 



Xh „ 



scattering suggested that no kinematical shift was neces- 
sary, i.e., Xh = x. This assumption has been recently 
questioned by a new analysis |2j|. Specifically, Tung et 
al. note that the heavy quarks in the hadron come pre- 
dominantly from gluons splitting into quark-antiquark 
pairs. Hence the heavy quark ft initiating the hard pro- 
cess must be accompanied by the unobserved h in the 
beam remnant. When both ft and ft are present, the 
hadron's light-cone momentum fraction carried by the in- 
coming parton cannot be smaller than x (l + 4Af 2 /Q 2 ), 
which is larger than the minimal momentum fraction 
imin = x allowed by the single-particle inclusive kine- 
matics. The factor of 4M 2 arises from the threshold con- 
dition for ft and ft. This effect can be accounted for by 
evaluating the flavor-excitation cross sections at the scal- 
ing variable Xh = % (l + 4M 2 /Q 2 ) . 

In brief, the rule proposed in Ref. |^| is to use Xa = 
x (l + 4M 2 /Q 2 ) in flavor-excitation processes (Fig. ||a) 
and Xa = x in flavor-creation processes (Figs. ||b and 
||c) when calculating inclusive cross sections. However, 
to correctly describe the differential distributions of the 
final-state hadron, we have to generalize the above rule 
for semi-inclusive observables. This generalization is dis- 
cussed in Appendix |C], where the proper scaling variable 
for fully differential finite-order cross section s is f ound to 
be Xh = x (1 + M 2 / (z(l - z)Q 2 )) (cf. Eq. (IcTl]) ). 



It is important to note that this prescription is not an 
approximation; it correctly accounts for the full mass de- 
pendence ||l7). It also tremendously reduces the com- 
plexity of flavor-excitation structure functions, as they 
are given by the light-quark result. In the specific case 
considered here, the heavy quark mass in the S-ACOT 
scheme should be retained only in the 7* + G — > ft + ft 
subprocess (Fig. ||b). Another important consequence 
will be discussed in Section 0, where we show that the 
S-ACOT scheme leads to a simpler generalization of the 
gT-resummation to the mass-dependent case. 



C. The scaling variable 

Finally, we address the issue of the most appropriate 
variables Xa (a = G,u,d,s, ...) in the convolution inte- 
gral (Q). In a massless calculation, Xa are just equal 
to Bjorken x, since all momentum fractions £ between 
x and unity are allowed by energy conservation. This 
simple relation does not hold in the massive case. For 
instance, in the charged-current heavy quark production 
W + q — > h, where h is present in the final, but not 
the initial, state, a simple kinematical argument leads 
to the conclusion that the longitudinal variable in the 
flavor-excitation processes should be rescaled by a mass- 
dependent factor, as Xh = x(l + M 2 /Q 2 ) [ |49[ . 

In the flavor-excitation subprocesses of the neutral- 
current heavy quark production (e.g., 7* + h — > ft), typi- 
cally no rescaling correction was made. The presence of a 
heavy quark in both the initial and final states of the hard 



III. MASSLESS TRANSVERSE MOMENTUM 
RESUMMATION 

We now turn to the differential distributions of the 
heavy-flavor cross sections. Specifically, we consider the 
production of a heavy-quark hadron H via the process 
e(£) + A(p A ) -> H(p H ) + e(£') + X. This reaction is il- 
lustrated in Fig. |] for the specific case when A is a pro- 
ton. In much of the discussion, we will find it convenient 
to amputate the external lepton legs and work with the 
photon-hadron process 7*(g) + A(pa) — * H{ph) + A" in 
the photon-hadron cm. frame. Given the conventional 
DIS variables Q 2 = —q 2 and x = Q 2 /(2pA ■ q), as well 
as the Lorentz invariant S e a = (£ + pa) 2 , we decompose 
the electron-level cross section into a sum over the func- 
tions A p (ip, tp) of the lepton azimuthal angle tp and boost 
parameter ip = cosh -1 (2xS e AQ~ 2 — l) 0, |5l|| : 



da{e + A^e + H + X) 
dxdQ 2 dpn 



00 ^>2 v p(q>Pa,Ph) 

p 

x A p (ip,ip). 



(4) 



This procedure is nothing else but the decomposition 
over the virtual photon's helicities [Hi); hence 

it is completely analogous to the tensor decomposition 
familiar from the inclusive DIS. As a result of this proce- 
dure, the dependence on the kinematics of the final-state 
lepton is factorized into the functions A p (tp,(p), while 
the hadronic dynamics affects only the functions V p . In 
parity-conserving SIDIS, the only contributing angular 



functions are 



where 



Ai(ip,ip) = 1 + cosh 2 ip, 

A 2 {^,ip) = -2, 

A 3 (tp,ip) = — cosipsinh2ip 7 

A4,(ip,<p) = cos2(/?sinh 2 if). (5) 

In Section |l| we found that the ACOT prescrip- 
tion resums logarithms of the form ln(M 2 /Q 2 ). For 
the inclusive observables, this procedure provides ac- 
curate predictions throughout the full range of x and 
Q 2 . More differential observables may contain additional 
large logarithms in the high-energy limit. In particu- 
lar, we already mentioned the logarithms of the type 
(q T 2 )agln m (gf7Q 2 ), < m < 2n - 1, which appear 
when the polar angle Oh of the heavy hadron H in the 
7* A cm. frame becomes small (cf. Fig. [j]). Here we 
chose the z-axis to be directed along the momentum q of 
the virtual photon 7*. When M 2 <§; Q 2 , the scale qr is 
related to Oh as 



hence 



q 2 T = Q 2 



lim q T = Q 



- 1 



1 — cos Oh 
1 + cos Oh 



1 



) 2 

7 H 



0. 



(6) 



(7) 



The resummation of these logarithms of soft and 
collinear origin can be realized in the formalism by 
Collins, Soper, and Sterman (CSS) {3J, || ff], H- The 
result can be expressed as a factorization theorem, which 
states that in the limit Q 2 ^> q 2 ^, {m 2 }, Aq CD the cross 
section is 



da(e + A^e + H + X) 



<z|«Q 2 



2S e A 



Ai{il>,<p) 



d 2 h 

Q 



e^ h WHA(b,Q lXl z) 

A-QCD 

Q 



(8) 



In this equation, 6 is the impact parameter (conjugate 
to qr), rf$ = dxdQ 2 dzdq^d(p, z = (pa ■ Ph)/(pa ■ g), 



and (To and Fi are constant factors given in Eq. (|B2j). 



As before, {m q } collectively denotes all quark masses, 
{m q } = m u ,md, M. At large Q 2 , the 6-space inte- 
gral in Eq. (||) is dominated by contributions from the 
region b 2 < l/Q 2 . In this region, the hadronic form fac- 
tor WnA(b, Q, x, z) can be factorized in a combination 
of parton distribution functions f a /A(£,a, Hf), fragmenta- 
tion functions Dnibi^b, Mf), and the partonic form factor 



W 



bit ■ 



W H a (b,Q,x,z) 



E 

a .. b 



1 d^a 
Ha 



1 d& 

6 



a 5 f^F J: 



,X,Z,fX F ) 



(9) 



W ba (b, Q, x, z, n F ) = J2 e 2 e ~s(b,Q, Cl ,c 2 ) 



z,ix F b;^-\CJ) a ( 



C 2 



(10) 



Here x = x/^ a , z = The indices a, b in Eq. (||) are 

summed over all quark flavors and gluons; the summation 
over j in Eq. (JITj) is over the quarks only. The fractional 
charge of a quark j is denoted as e 2 e 2 . The parton distri- 
butions and fragmentation functions are separated from 

the partonic form factor Wba at the factorization scale 
fiF- The Sudakov factor S(b, Q, C\, C2) is an all-order 
sum of logarithms ln m {q^/Q 2 ). It is given by an inte- 
gral between scales C\jb 2 and C 2 Q 2 (where G\ and C2 
are constants of order 1) of two functions A (as(ft); C\) 
and B{as{p)\ Ci, C2) appearing in the solution of equa- 
tions for renormalization- and gauge-group invariance: 



S 



C- 2 Q- 



dfL 2 



Cf/b 2 M 



111 



C 2 2 Q' 



^(a s (jS);Ci) 



+B{a s {jx);C 1 ,C 2 ) 



(11) 



The functions C m , C out contain perturbative correc- 
tions to contributions from the incoming and outgoing 
hadronic jets, respectively. To evaluate the Fourier- 
Bessel transform integral, WnA(b,Q,x,z) should be 
also defined at b > 1 GeV -1 , where the perturba- 
tive methods are not trustworthy. The continuation of 
Wha (b, Q, x, z) to the large-6 region is realized with the 
help of some phenomenological model, as discussed, e.g., 
in Refs. §|, || || . 

As noted above, the resummed cross section in Eq. (||), 
which we shall label as er~ , is derived in the limit q\ -C 
Q 2 . In the region q\ > Q 2 , the standard finite-order 
(FO) perturbative result, <jfo, is appropriate. While 
and ijfo represent the correct limiting behavior, we can- 
not simply add these two terms to obtain the total cross 
section, (Jtot, as we would be "double-counting" the con- 
tributions common to both terms. 

The solution is to subtract the overlapping contribu- 
tions between cr^ and ofo- This overlapping contribu- 
tion (the asymptotic piece a asy) is obtained by expand- 
ing the 6-space integral in out to the finite order of 
o~fo- Thus, the complete result is given by 



da 



TOT 



do 



doFO da asy 



d$> 



d4> 



d<f> 



(12) 



At small qr, where terms ln m (g|,/Q 2 ) are large, opo 
cancels well with <jasy, so that the total cross section is 
approximated well by the 6-space integral: otot ~ o 



At q\ > Q 2 , where the logarithms are no longer domi- 
nant, the 6-space integral cancels with a asy, so that 
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VTOT ~ + (T F &ASY 




(b) 

Figure 3: Balance of various terms in the total resummed cross section doTOT / dqr: (a) away from the threshold (Q ^> M); 
(b) near the threshold (Q ~ M). In each plot the thick curves correspond to the "active" cross section (TOT, FO, W or ASY), 
and the thin curves correspond to the other three cross sections. 



the total cross section is equal to crpo up to higher order 
corrections: <jtot ~ &fo- This interplay of c^, ctfo, 
and a asy in (Jtot is illustrated in Fig. ||a. 

As we will be referring to these different terms fre- 
quently throughout the rest of the paper, let us present 
a recap of their roles. 

• cr^ is the small-gr resummed term as given by 
the CSS formalism in Eq. (||); sometimes called 
"the CSS term" [Q. This expression contains 
the all-order sum of large logarithms of the form 
ln m (q^/Q ), which is presented as a Fourier-Bessel 
transform of the 6-space form factor W(b, Q, x, z). 
It is a good approximation in the region < Q 2 . 

• ofo is the finite-order (FO) term; sometimes called 
"the perturbative term". It contains the complete 
perturbative expression computed to the relevant 
order of the calculation n. As such, this term con- 
tains logarithms of the form In™ (q^/Q 2 ) only out 
to m = In — 1. It also contains terms that are not 
important in the limit q\IQ 2 — > 0, but dominate 
when q\ ~ Q 2 . Hence, it provides a good approxi- 
mation in the region q^^i Q 2 - 

• &asy is the asymptotic (ASY) term. It contains the 
expansion of er~ out to the same order n as in apo- 
As such, this term contains logarithms of the form 
ln m (q^/Q 2 ) only out to m = 2n— 1. It is precisely 



what is needed to eliminate the "double-counting" 
between the a^, and <tfo terms in Eq. (p^). 

• <Jtot is the total (TOT) resummed cross section; 
sometimes called "the resummed term". It is con- 
structed as o~tot = + &fo — &ASY- In the 
region q\ <C Q 2 , a asy precisely cancels the large 
terms present in the <jfo contribution, so that 
otot ~ cr~. In the region q\ > Q 2 , a asy ap- 
proximately cancels the term leaving ofo as 
the dominant representation of the total cross sec- 
tion: otot ~ &fo- Hence, when calculated to a 
sufficiently high order of as, &tot serves as a good 
approximation at all gr- 
in a practical calculation in low orders of PQCD, one 
may want to further improve the cancellation between 
cr^ and a asy at q\ > Q 2 . This improvement can 
be achieved by introducing a kinematical correction in 
these terms that accounts for the reduction of the al- 
lowed phase space for the longitudinal variables x and z 
at non-zero qr- The purpose of this kinematical correc- 
tion is quite similar to the purpos e of t he inclusive scaling 
variable discussed in Subsection [lIC| : it removes contri- 
butions from the unphysically small x and z that make 
the difference er^> — a asy non-negligible as compared to 
(7 fo . Note that the resummed cross sections with and 
without the kinematical correction are formally equiva- 
lent to one another up to higher-order corrections. Fur- 
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ther discussion of this issue can be found in Appendix |C], 
which introduces the kinematical correction for the re- 
summed heavy-quark qr distributions. 

IV. EXTENSION OF THE CSS FORMALISM TO 
HEAVY-QUARK PRODUCTION 

In the previous Section, we presented a procedure for 
the resummation of distributions dajdq\ in the limit 
when Q 2 is much larger than all other momentum scales, 
Q 2 3* Qti { m q}- We now are ready to discuss its ex- 
tension to the case when the heavy-quark mass is not 



negligible. For simplicity, we again assume that only 
one heavy flavor h has the mass M comparable with Q: 
Q 2 ~ M 2 3> Aq C£) . The generalization for several heavy 
flavors can be realized through the conventional sequence 
of factorization subschemes, in which the heavy quarks 
become active partons at energy scales above their mass, 
and are treated as non-partonic particles at energy scales 
below their mass. 

We start by rewriting Eq. (Q) in a form analogous to 
Eq. (4.3) of Ref. |3(|, where the form factor W was given 
for the Drell-Yan process: 



~ , . ^ . ^ — ^ 9 ~Pr5i ou ^ (if -i C\ I ~^^ n (if -i C 1 

W H a [b, Q, x,z) = e J 3* H/j \z, b, {m g }, — I 3» j/A I x, b, {m q }, — 



x exp<, 



dp 2 



[ ^ )A ( a s (A); \ — h Ci ) + B ( a s (p,); I JL } ; C x , C 2 



(13) 



r 



Here the function ^ , A (x, b, {m g } , Ci /C2) describes 
contributions associated with the incoming hadronic jet. 

As illustrated in Appendix [A|, S^j/A i x i b, {wiq}, C1/C2) 
is related to the fcx-dependent parton distribu- 
tion &J/a ( x i ^t, {mq})- Similarly, the function 

^"hh (2, b, {m q }, C1/C2) describes contributions asso- 
ciated with the outgoing hadronic jet |35|]. It is 
related to the ^-dependent fragmentation function 
^"h% ( z 1 , { m q}) ■ The functions A and B are the same 
as in Eq. (|ll]), except that now they retain the depen- 
dence on the quark masses {m q } = m u , md, m s , ...,M. 

Eq. (j|) presents a special case of Eq. ( |l3|) . It is valid 
at short distances, i.e., when 1/6 is much lar ger than any 
of the quark masses m q . In contrast, Eq. (113) is valid 



at all b 



As shown in Ref. [|45fl, the transition from 



Eq. (|13|) to Eq. fe|) is possible because the functions @* 



j/A 



and ^H/j factorize when b\jb 2 ^> {? 



Cx 
C 2 



& 'j/A [x,b,{m q }, 



E 



6> 



—r-OUt 



fa/ A 



m q 



• b,{ m q},7T 



6 ^ & 



ff/6 



m„ 



pout 
L b/j 



Z, /ip6: 



Co 



(14) 



Here we introduced a frequently encountered constant 
b = 2e,-^ E w 1.123. We see that the form-factor Wha 
is well-defined both for non-zero quark masses and in 



the massless limit. Hence, it does not contain negative 
powers of the quark masses or logarithms In (m q /Q), with 
the exception of the collinear logarithms resummed in the 
parton distributions and fragmentation functions. 

We will now argue that the factorization rule similar 
to Eq. (H) should also apply in heavy-flavor production 
when M^is not negligible compared to b 2 /b 2 . Indeed, 

ifi out 

the factorization of the functions j/ A an d ^H/j m the 
limit b\jb 2 3> {m 2 } JiB]] closely resembles the factoriza- 
tion of the inclusive DIS structure functions in the limit 
Q 2 > {m 2 q } §§, IgUp, H, H. In both cases the fac- 
torization occurs because the dominant contributions to 
the cross section come from "ladder" cut diagrams with 
subgraphs containing lines of drastically different virtu- 
alities. More precisely, the leading regions in such di- 
agrams can be decomposed into hard subgraphs, which 
contain highly off-shell parton lines; and quasi-collinear 
subgraphs, which contain lines with much lower virtual- 
ities and momenta approximately collinear to p A (in the 

case of F h / A (x, Q 2 ) or @* j / A ) or (in the case of £? H /j). 

271 Out 

In the functions & ji A and ^jj/jj additional soft gluon 
subgraphs are present, but they eventually do not affect 
the proof of the factorization []45| . The hard subgraphs 
contribute to the inclusive coefficient function C/j/ a in 
Eq. (§), as well as functions Cf, or C°J* in Eq. @. 
The quasi-collinear subgraphs, which are connected to 
the hard subgraphs through one on-shell parton on each 
side of the momentum cut, contribute to the PDF's (in 
the inclusive DIS and SIDIS) or FF's (in SIDIS). 

The hard subgraphs are characterized by typical trans- 
verse momenta k\ > fi F ^> Aq CD , while the PDF's and 
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FF's are characterized by transverse momenta kp < p? F . 
The factorization scale [ip is of order Q in the inclu- 
sive DIS structure functions and bo /b in the functions 

^J/A and ^°H/j- As discussed in Section ||, the fac- 
torization in the inclusive DIS can be extended to the 
case when Q is comparable to the heavy-flavor mass M , 
Q 2 ~ M 2 S> Ag C£) . Given the close analogy between 
the inclusive DIS structure functions and the functions 

— ~zn — r^OUt 

™j/Ai ^H/j ■> !t ls natural to assume that the latter 
factorize when b\/b 2 ~ M 2 > Aq C£ > as wen: 
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(15) 



The main difference between Eqs. ([L4|) and ( p. 5|) is con- 
tained in the functions Cj/ Q and C£J- , which now explic- 
itly depend on M. These functions can be calculated 
according to their definitions given in Ref. [Q. The 
unrenormalized expressions for the C-functions contain 
ultraviolet singularities. To cancel these singularities, 
we introduce counterterms according to the procedure 
described in Section [[J that is, graphs with internal 
heavy-quark lines are renormalized in the MS scheme 
if fip ~ bo/b > M and by zero-momentum subtraction if 
bo/b < M. This choice leads to the explicit decoupling 
of diagrams with heavy quark lines at b > bo/M. In 
particular, the decoupling implies that contributions to 
Eq. (|T^ ) with j, a, or b equal to h are power-suppressed 
at b > b /M. 

We now consider other sources of the dependence on M 
in da/d<&. Firstly, according to Eq. (|lj), there is a depen- 
dence on M in the Sudakov functions A(as(p); p>/M; C\) 
and B(ag(ji);p,/M;Ci,C2). Due to the decoupling, the 
mass-dependent terms in the Sudakov factor vanish at 
b > bo/M, except for perhaps terms of truly nonpertur- 
bative nature, like the intrinsic heavy quark component 
|p7| . As mentioned above, in this paper such nonpertur- 
bative component is ignored. Secondly, there may also 
be mass-dependent terms in the finite-order cross sec- 
tion, which are not associated with the leading contri- 
butions resummed in the W-term: those are the terms 
that contribute to the remainder in Eq. (||). The terms 
of both types are correctly included in do~pop / d<&. In- 
deed, the terms of the first type appear in all three terms 



da^/d^, dapo/d&, and da^sY /d&. Two out of three 
contributions (in do^/d<& and d&ASY / d&, or dapo/d^ 
and doAsy/d^) cancel with one another, leaving the 
third contribution uncancelled in daror/d^. The terms 
of the second type are contained only in dapo/d<&, so 
that they are automatically included in dapor / 'd$> . 

The treatment of the massive terms simplifies more 
if we adapt the S-ACOT factorization scheme, in which 
the heavy quark mass is set to zero in the hard parts 
of the flavor-excitation subprocesses. As a result, M 
is neglected in the flavor-excitation contributions to the 
hard cross section opo, asymptotic term ctasy, and C- 
functions in the W-term. The mass-dependent terms are 
further omitted in the perturbative Sudakov factor S. At 
the same time, all mass-dependent terms are kept in upo, 
o~asy, and C-functions for gluon- initiated subprocesses. 

As we will demonstrate in the next section, in this 
prescription the cross section (Jtot resums the soft and 
collinear logarithms, when these logarithms are large, 
and reduces to the finite-order cross section, when these 
logarithms are small. In particular, at Q ~ M the finite- 
order flavor-creation terms approximate well the heavy- 
quark cross section. Hence we expect that upor repro- 
duces the finite-order flavor-creation part at Q ~ M 
(Fig. ||b). For this to happen, the flavor-excitation 
cross section should cancel well with the subtraction 
cx ln(/ii?/M) from the flavor-creation cross section; and 
cr^ should cancel well with (T^sy.We find that these can- 
cellations indeed occur in the numerical calculation, so 
that at Q « M apor agrees well with the flavor-creation 
contribution to o~po- Similarly, cttot reproduces the 
massless resummed cross section when Q M (Fig. ||a). 
It also smoothly interpolates between the two regions of 
Q. 

To summarize our method, the total resummed cross 
section in the presence of heavy quarks is calculated as 



(7(7 



TOT 



da 



da 



FO 



(/<!> 



(/<!> 



d(JAS\ 

d<5> 



(16) 



i.e., using the same combination of the W-term, finite- 
order cross section, and asymptotic cross section as in the 
massless case. All three terms on the r.h.s. of Eq. ( |l^ ) 
are calculated in the S-ACOT scheme. The W-texm is 
calculated as 



da(e + A-> e + H + X)\ _ <tqF\ Aj^tp) 

S e A 2 



(/<!> 



d 2 h 

(2tt) 2 



e^ Tb W HA {b,Q,M,x,z), 



(17) 



where the form-factor WhaQ), Q, M, x, z) is 
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and 



r c * Q an 2 

S ba (b,Q,M) = / 4r 



A(a s (-p); C,) In ( ^ 



+B(a s (7Z);Ci,C 2 ) 



+ ^ P (6,Q,M). (19) 



As in the factorization of inclusive DIS structure func- 
tions (cf. Section |J), we find it useful to replace Bjorken 
x by scaling variables 



Xh=x[l + 



M 



f2 



z(l-z)Q 2 , 

in (Tfo for the flavor-excitation subprocesses, and 

M 2 + z 2 q\" 
z(l-z)Q 2 



Xh = x ( 1 + 



(20) 



(21) 



in o and oasy • The purpose of these scaling variables 
is to enforce the correct threshold behavior of terms with 



incoming heavy quarks. Eqs. (|20j) and ( |2l| ) are derived 
in detail in Appendix ^|. 



V. MASSIVE RESUMMATION FOR 
PHOTON-GLUON FUSION 

We now analyze contributions to the total resummed 
cross section daror/d^ from the 0(oP s ) heavy-flavor ex- 
citation subprocess 7*(g) + h(p a ) — > h(pb) (Fig. ||a) and 
(D(as) photon-gluon fusion subprocess 7*(g) + G(p a ) — > 
h(p b ) + h(p s ) (Fig. |b). Since we work in the S-ACOT 
scheme, only the 0(as) fusion subprocess retains the 
heavy quark mass, so that we concentrate on that pro- 
cess first. The expression for the 7*/i contribution, which 
is the same as in the massless case, is given in Eq. (Bl). 
In the following we outline the main results, while details 
are relegated to Appendices. 



"recycle" as much of the results from Refs. |pl| , p2\ |5l|| as 
possible. 

Throughout the derivation, the mass of the incident 
hadron will be neglected: p\ = 0. We will use the stan- 
dard DIS variables x, Q 2 , and z, defined by 



2pA ■ q 



2_ 2. _Pa-Ph 



PA ■ q 



Q=-q 



(22) 



Since we will be interested in the transverse momentum 
distributions (or equivalently, the angular distributions), 
we next define the transverse momentum in a frame- 
invariant manner. The four-vector gf of the transverse 
momentum must be orthogonal to both of the hadrons, 
so that we have the conditions qfPA = and qt -pn = 0. 
In the massless case, gf is simply defined by subtracting 
off the projections of the photon's momentum g M onto pa 
and ph- This is slightly modified in the massive case to 
become 



It 



q' 



PH ■ q 



- M 



Pa ■ q 



PA-PH H {PA-Ph) 2 



u 

Pa 



PA-q ft 

Ph- 

Pa - Ph 



(23) 



Here Mr denotes the mass of the heavy hadron. We find 
for q\ = 



-qfqtp. 



2 

q T 



Q 



2 PH ■ q Mjj 

z z 2 



(24) 



The kinematical variables at the parton level can be in- 
troduced in an analogous manner. Let £ a denote the frac- 
tion of the large '— ' component of the incoming hadron's 
momentum pa carried by the initial-state parton a (i.e., 
£a = Pa Ip~a)'? an d 6> denote the fraction of the large 
'+' component of the final-state parton's momentum p 
carried by the outgoing hadron H (i.e., = Pn/pt)- 
We also assume that relates the transverse momenta 
of b and H, as (j>t)h = £b(j>T)b. Since all incoming par- 
tons are massless in the S-ACOT factorization scheme, 
we find the following relations between the hadron-level 
variables x, z, qr and their parton-level analogs x, z, qr- 



(25) 



2 (j>a ■ q) & ' 



Mass-Generalized Kinematical Variables 



Our approach will be to first generalize the kinematical 
variables from the massless resummation formalism to 



3 We remind the reader that the analysis is performed in the 7* A 
cm. frame, where the incident hadron moves in the — z direction. 
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It 



= (Pa ■ Pb) 

(Pa ■ q) 



(26) 
(27) 



where in the derivation of Eq. (|27|) we used the first equal- 
ity in Eq. @. 

If we introduce a massive extension of q\ called q\ and 
defined by 



q T = q T 



M 2 



(28) 



then the form of q\ in terms of the Lorentz invariants is 
identical to the massless case: 



f T = Q 2 + 2 



Ph ■ q 



z 



(29) 



We also generalize the usual Mandelstam variables 
{§, t, u} to what we label the "mass-dependent" Mandel- 
stam variables {J, ti,ui}: 

s = (q+Paf, (30) 
h = t- M 2 = (q- p h f - M 2 , (31) 
ux = u- M 2 = ( Pa - Ph ) 2 - M 2 . (32) 

By using the variables q^,s,ti, and u\ instead of their 
counterparts q T , s, t, and u, we shall be able to cast 
many of the massive relations in the form of the mass- 
less ones. For example, the expressions for the "mass- 
dependent" Mandelstam variables {s, in terms of 
the DIS variables can be written as 



141 



-Q 



Q 2 (z-\)-q%z = -< 



,(1-2-) 



(33) 
(34) 
(35) 



Note how we made use of the generalized transverse mo- 
mentum variable q\. These relationships have the same 
form as their massless counterparts. As a result, the 
denominators of the mass-dependent propagators, which 
are formed from the invariants s, t\, and ui, retain the 
same form as the denominators of the massless propaga- 
tors, which are formed from the invariants s", t, and u. 



B. Relations between {Eh, cos 6n}in the 7*A 
cm. frame and {z,q^} 

It is useful to convert between the final-state energy 
Eh, polar angle 8h and the Lorentz invariants {z,q^}. 
Given the j*A cm. energy W 2 = (q+pA) 2 = Q 2 (l — x)/x 

one easily finds the follow- 



andp= \p H \ = \/E% 
ing constraints on Eh, p, and cos 6h 



M 2 



M H <E H < 



ir- 



< v < 



(l-*& 



(36) 
(37) 



q T /W 



10' 



10" 



10 1 



10" 



10" 



Current region 



Target region 




25 50 75 100 125 150 175 



Figure 4: Plots of qr/W vs. Oh at various values of A = 
Mh/Eh — 0.999 (lower curve), A = 0.5 (middle curve), and 
A = 0.001 (upper curve). 



and 

-I<cos0fl-<1. (38) 
Given Eh and cos Oh , we can determine z and q\ as 



1 

W 



(E H +p cos Oh) 



q% 



(P TJ 



w 



p 2 (1 - cos 2 H ) 
(E h +p cos H ) 



From Eqs. 



8h the bounds on z can be found as 



W 2 



< z < 1. 



(39) 
(40) 

(41) 



Note that the first equality in Eq. ( fl0| ) identifies qr as 
the the transverse momentum of H rescaled by the final- 
state fragmentation variable z. Hence qr can be also 
interpreted as the leading-order transverse momentum 
of the fragmenting parton. Similarly, qr = Mt/z can 
be interpreted as the rescaled transverse mass Mt of the 
heavy quark. It also follows from Eqs. ([?9 40) that the 
two-variable distribution with respect to the variables z 
and qr coincides with the two-variable distribution with 
respect to Eh and Oh'- 



da 



da 



dxdQ 2 dzdqT dxdQ 2 dEnd0n 



(42) 



As a result, the distributions in the theoretical variables 
z and qx are directly related to the distributions in Eh 
and Oh measured in the experiment. 

Despite the simplicity of the relation (|4^), z and qr 
are quite complicated functions of Eh and cos Oh indi- 
vidually. This feature is different from the massless case, 
where there exists a one-to-one correspondence between 
qr and cos Oh for the fixed 7* A cm. energy W: 



W 



q T 



\Mh=0 



W 2 



Qt 



(43) 
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This relationship does not hold in the massive case, in 
which one value of qx corresponds to two values of cos Oh- 
Indeed, Eq. (En) can be expressed as 



a% (1 - A 2 ) (1 - cos 2 H ) 



W 2 



(i + vr~^coso H y 



(44) 



where, according to Eq. (|3q) , the variable A = AIh/Eh 
varies in the following range: 



2M H 



W (1 + M H /W 2 ) 
Eq. (EJ) can be solved for cos H as 

1 



< A < 1. 



(45) 



cos Vh 



(q 2 + W 2 ) v / T^A 2 ~ 



When the energy Eh is much larger than Mh (A — * 0) 
the solution with the "+" sign in Eq. ( ^6| ) turns into the 
massless solution ([I||). The solution with the "— " sign 
reduces to cos Oh = — 1. 

The physical meaning of the relationship between qr 
and cos Oh can be understood by considering plots of 
qr/W vs. Oh for various values of A (Fig. ||). Let us 
identify the current fragmentation region as that where 
cosOh is close to +1 (Oh = 0), and the target fragmenta- 
tion region as that where cos Oh is close to —1 (Oh = 7r). 
Firstly, qx = if cos Oh — 1 or cos Oh = — 1 . Secondly, 
near the threshold (A — ► 1) the ratio qr/W is vanishingly 
small and symmetric with respect to the replacement of 
Oh by (it — Oh)- Thirdly, as A decreases, the distribution 
qr/W vs. Oh develops a peak near Oh = 180°. This 
peak is positioned at cos Oh = — vl — A 2 , and its height 

is qr/W = (l - A 2 ) 1/2 /A 2 . For H < 180°, the distribu- 
tion rapidly becomes insensitive to A; more so for smaller 
Oh- 

In the limit A — > 0, the peak at Oh = 180° turns 
into a singularity. This singularity resides at the point 
z = and corresponds to hard diffractive hadroproduc- 
tion. The analysis of this region requires diffractive par- 
ton distribution functions |64|, ||[ j6(| Q and will not 
be considered here. For Oh ^ 180°, one recovers a one- 
to-one correspondence between qr/W and cos Oh of the 
massless case. We see that there is a natural relationship 
between and cos Oh, which becomes especially simple 
in the massless limit. In the following, we concentrate on 
the limit qr — > and z ^ 0, which corresponds to the 
current fragmentation region 6h — * 0. 



C. Factorized cross sections 

Next, we consider the factorization of the hadronic 
cross section. Given the hadron-level phase space ele- 
ment d$ = dxdQ 2 dzdq T dip and its parton-level analog 



d$ = dxdQ 2 dzdq^df>, all three terms on the r.h.s. of 
Eq. (JLq) can be written as 



E 



da 



x D H / b 6 



1 d£ b f 1 d£ a 



(Ip 



fa/ A I £a, 



(i F 



dd fxa z q T (i F M 



-q^±WJ(l-\ 2 )(q 2 T + W 2 )-q 2 T \. (46) ment \M\' 



d$ Q' Q 1 Q 1 

Let us first consider the finite-order cross section 
dapo/d^ . The explicit expression for this cross sec- 
tion at the lepton level can be found in Appendix [B|. We 
are interested in extracting the leading contribution in 
this cross section in the limit Q — > oo with other scales 
fixed. Specifically, we concentrate on the behavior of the 
phase-space 5— function that multiplies the matrix ele- 



^ oc 6 (s + t + u + Q 2 - 2AI 2 ) \M\ 2 
= S (s + h +ui + Q 2 ) \M\ 2 



S \ \ - 

x 



1 

-- 1 

z 



f,)\M\. (48) 

Here we used the mass-generalized variable q\ introduced 
in Eq. (|2^) . Note that in terms of the variables x, z, and 
q 2 ^ this expression takes the same form as its massless 
version. In the limit Q — > oo, and x,~z, and qr fixed, the 
5-function can be transformed using the relationship 



lim S (yiy2 - e) 

£^0 



This transformation yields 
lim 6 (s + t + 



S(yi) S(y 2 ) 
[V2\+ [yi}+ 
log(e)6(y 1 )S(y 2 ) 



(49) 



Q 2 - 2M 2 ) 



6(1 



6(1 



log 



3t 



6(1 -x) 6(1 - z). 



(50) 
(51) 



This asymptotic expression for the (5-function is exactly 
of the same form as in the massless case up to the re- 
placement q\ — > q\. 

Furthermore, in the above limit the matrix element 
\M\ 2 itself contains singularities when Q 2 3> q\. In 
particular, the largest structure function V\ in the 7*G- 
fusion subprocess (cf. Eq. flB4D) contains contributions 
proportional to 



1 



1 



(M 2 -t)(M 2 -u) hm 



1 



and 



M 2 M 2 
- oc — 
It 



z^M 2 



t{u\ 



(Pg 2 + M 2 ) 



2 ■ 



(52) 



(53) 
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When M is not negligible, these contributions are finite 
and comparable with other terms. However, in the limit 
when both M and qr are much less than Q, the terms of 
the first type diverge as The terms of the second 

type vanish at qr 7^ and yield a finite contribution at 
qx = 0. These non-vanishing contributions are precisely 
the ones that are resummed in the VF-term; in the total 
resummed cross section axor, they have to be subtracted 
in the form of the asymptotic cross section a asy to avoid 
the double-counting between <jfo and a^. 

To precisely identify these terms, we calculate them 
from their definitions, as described in Appendix [a|. Since 
the 0{as) 7*G subprocess is finite in the soft limit, it 
contributes only to the function CJ^, G (x, /j,Fb,bM) and 
not to the Sudakov factor. The O(asfn) coefficient in 
this function is 

Qin(l) ( 



:^>(x,fi F b,bM) = T R x(l - x) (1 + Cl (bM)) 



Pt%(x) (c (6M)-ln(^) 



(54) 



if fiF > M, and 

C h/G bM ) 



C Z/G ^Fb, bM) 



fj, F >M 



fJ-F 

M 



(55) 



if fip < M. 



tion: P, 



(i) 



Here P, 

= T, 



(i) 



h/G^I ls ^ ne MS splitting func- 
h/G (0 =T R (l-2£ + 2e) , with T R = 1/2. The 
functions co(bM), and a(bM) denote the parts of the 
modified Bessel functions K (bM) and bM K\(bM) that 
vanis h wh en b -C 1/M. They are defined in Eqs. (A16) 
and (A19), respectively. 

We now have all terms necessary to calculate the 



combination 



/v>m(0) 

y^h/h 



® fh/A){x) 



(C/G ® fG/A)(x), 



which serves as the first approximation to the function 

ZP'h/A ( x > b, M, C1/C2). We find that this combination 
possesses two remarkable properties: it smoothly van- 
ishes at up = b\jb 2 <C M 2 and is differentiable with 
respect to ln(/ij?/M) at the point /xf = M. As a result, 
the form factor W{b,Q,x,z) for the combined 0{a s ) 
flavor-excitation and 0{a l s ) flavor-creation channels is 
a smooth function at all b, which is strongly suppressed 
at b 2 b^/M 2 . The physical consequence is that, for 
a sufficiently heavy quark, the 6-space integral can be 
performed over the large-6 region without introducing an 
additional suppression of the integrand by nonperturba- 
tive contributions. We use this feature in Section [vj, 
where we calculate the resummed cross section for bot- 
tom quark production, which does not depend on the 
nonperturbative Sudakov factor. 

Finally, by expanding the form-factor Wha in a series 
of as /it and calculating the Fourier-Bessel transform in- 
tegral in Eq. ( |l7| ) , we find the following asymptotic piece 
for the 7*G fusion channel: 



' da(e + G -> e + h + h) 



gqFi as 

ASY 4nSeA 7T 



xA 1 (tp,cp)6(l-z) 



P h/G@) , M 2 x{l-x) 



~2 

q T 



^A 

q T 



(56) 



When Q ~ M, da asy I dQ, which is a regular function at 
all qr, cancels well with da~/d§. In the limit Q — > 00, 
da asy /d<& precisely cancels the asymptotic terms that 
appear in the finite-order cross section doFo/d'b. 

VI. NUMERICAL RESULTS 

In this Section, we apply the resummation formalism 
to the production of bottom quarks at HERA. The calcu- 
lation is done for the electron-proton cm. energy of 300 
GeV and bottom quark mass M = 4.5 GeV. For simplic- 
ity we assume that the masses of the B-hadrons coincide 
with the mass of the bottom quark M. We also neglect 
the mixing of photons with Z°-bosons at large Q. 

In the following, we discuss polar angle distributions 
in the j*p frame for x = 0.05 and various values of Q. 
The cross section is calculated in the lowest-order ap- 
proximation as discussed in Section [y]. 4 The calculation 
was realized using the CTEQ5HQ PDF's [§g| and Pe- 
terson fragmentation functions |7(J with e = 0.0033 [||. 
The finite-order cross section daFo/dQ and asymptotic 
cross section da asy /d& were calculated at the factoriza- 
tion scale [If = Q- The scale- related constants in the 
W-term were chosen to be C\ — 2e~~ lE = bo and C2 = 1, 
and the factorization scale was [If = bo/b. The W-term 
included the 0{a s ) C-functions (x, /if&, C1/C2), 

(z, (Xfo, C1/C2) and 0(a s ) initial-state function 

C™,Q (x, (iFb, bM). In addition, it included the pertur- 
bative Sudakov factor (|Tl|), unless stated otherwise. The 
Sudakov factor was evaluated at order 0(as), which was 
sufficient for this calculation given the order of other 
terms. The functions in the Sudakov factor were eval- 
uated as 



A( f i;C 1 ) = C J 



and 



3C F as{n) 



(57) 



(58) 



According to the discussion in Section 0, our calcu- 
lation ignores unknown nonperturbative contributions in 



4 The generalization of our approach to higher orders is straight- 
forward. The next-order calculation should include the 0(as) 
flavor-excitation and 0(a 2 s ) flavor-creation channels, which 
should appear together to ensure the smoothness of the form 
factor W(b, Q, M, x, z) and its suppression at b > 1/M. 
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Figure 5: The angular distributions of the bottom hadrons in the 7*p cm. frame at (a) Q = 5 GeV, (b) Q = 15 GeV, (c) 
Q = 50 GeV without the Sudakov factor, and (d) Q = 50 GeV with the Sudakov factor. At Q = 50 GeV, an additional cut 
Eh > 0.1(W/2) is made to suppress contributions at z < 0.1, i.e., from the region where the conventional factorization may be 
inapplicable. The plots show the finite-order cross section ofo (long-dashed line), the 6-space integral a~ (dot-dashed line), 
the asymptotic piece a asy (dotted line), and the full resummed cross section otot (solid line). 



the PF-term. In the numerical calculation, we also need 
to define the behavior of the light-quark PDF's at scales 
= ba/b < 1 GeV. Due to the strong suppression of 
the large-6 region by the M-dependent terms in the C- 
functions (cf. the discussion after Eq. ( A12 )), the exact 
procedure for the continuation of the PDF's to small u,p 
has a small numerical effect. We found it convenient to 
"freeze" the scale fip at a value of about 1 GeV by in- 
troducing the variable b* = bj yl + (b/b max ) 2 |3(J with 
bmax — ^0 GeV" 1 w 1.123 GeV" 1 . Other procedures 
|^6| , [5^] for continuation of WjiA(b, Q, x, z) to large val- 
ues of b may be used as well. Due to the small sensitivity 
of the resummed cross section to the region b 2 S> b 2 /M 2 , 
all these continuation procedures should yield essentially 
identical predictions. 

Fig. H demonstrates how various terms in Eq. jl6| ) 
are balanced in an actual numerical calculation. Near 
the threshold (Q = 5 GeV, Fig. ||a) the cross section 
d<JTOT / (dxdQ 2 d9 h ) should be well approximated by the 



0(as) flavor-creation diagram 7* + G — > h + h. We find 
that this is indeed the case, since the W^-term, which does 
not contain large logarithms, cancels well with its pertur- 
bative expansion duAsy/d^. As a result, the full cross 
section is practically indistinguishable from the finite- 
order term. 

At higher values of Q, we start seeing deviations from 
the finite-order result. Fig. ||b shows the differential 
distribution at Q = 15 GeV, i.e., approximately at 
Q 2 /M 2 = 10. At this energy, da T OT I '(dxdQ 2 dd H ) still 
agrees with dapo/ (dxdQ 2 d9H) at 9h ^ 10°, but is above 
do fo I (dxdQ 2 d6 h ) at 9h ^ 10°. The excess is due to the 
difference da ^ / '(dxdQ 2 d6 h) — da asy I '(dxdQ 2 d6 h) , i.e., 
due to the higher-order logarithms. 

Away from the threshold (Q = 50 GeV), 
daroT / (dxdQ 2 d6n) is substantially larger than the 
finite-order term at Oh ^ 10°, where it is dominated 
by da ^ / (dxdQ 2 dO h) ■ In this region, da fo I (dxdQ 2 dO h) 
is canceled well by da asy / (dxdQ 2 dO h) ■ Note, however, 
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Figure 6: Comparison of the massive and massless cross sec- 
tions at (a) Q = 15 GeV, and (b) Q = 50 GeV. The plots 
show the massive resummed cross section cttot (thick solid 
line); the massless resummed cross section apor (thin solid 
line); the massive finite-order cross section ctfo (thick dashed 
line); and the massless finite-order cross section ctfo (thin 
dot-dashed line). 



that contrary to the experience from the massless case, 
daFo/(dxdQ 2 d0H) and da asy / [dxdQ 2 dO h) are not sin- 
gular at 9h — * due to the regularizing effect of the 
heavy quark mass in the heavy-quark propagator at 
Oh ^ 3°. Figs. ||c and ^|d also compare the distribu- 
tions with and without the 0{a 1 s ) perturbative Sudakov 
factor, respectively. Note that at the threshold the flavor- 
excitation terms responsible for S are of a higher order 
as compared to the 0{pP s ) flavor-excitation and 0(a s ) 
flavor-creation terms. Correspondingly, near the thresh- 
old the impact of S is expected to be minimal. This 
expectation is supported by the numerical calculation, in 
which the difference between the curves with and with- 
out the 0(as) perturbative Sudakov factor is negligi- 
ble at Q = 5 GeV, and is less than a few percent and 
Q = 15 GeV. In contrast, at Q = 50 GeV the distribution 
with the 0(as) Sudakov factor is noticeably lower and 
broader than the distribution without it: at some val- 
ues of Oh, the difference in cross sections reaches 40%. 
The influence of the Sudakov factor on the integrated 



rate is mild: the inclusive cross section da /{dxdQ 2 ) cal- 
culated without and with the 0{as) Sudakov factor is 
equal to 330 and 320 fb/GeV 2 , respectively. Due to 
the enhancement at small Oh, these resummed inclu- 
sive cross sections are larger than the finite-order rate 
da FO /(dxdQ 2 ) « 260 fb/GeV 2 by about 25%. 

It is interesting to compare our calculation with the 
massless approximation for the j*G contribution. Fig. || 
shows the finite-order and resummed cross sections cal- 
culated in the massive and massless approaches. In con- 
trast to the massive cttot, the massless apor must in- 
clude the nonperturbative Sudakov factor S NP , which is 
not known a priori and is usually found by fitting to the 
data. To have some reference point, we plot the mass- 
less a TO T with S NP {b) = b 2 M 2 /b 2 r» 166 2 , so that, in 
analogy to the massive case, the region of b > b n /M in 
the massless W(b, Q, x, z) is suppressed. Since the heavy- 
quark mass has other effects on the shape of W(b, Q, x, z) 
besides the cutoff in the 6-space, we expect the shape of 
the massless and massive resummed curves be somewhat 
different. This feature is indeed supported by Fig. ||b, 
where at small Oh both resummed curves are of the same 
order of magnitude, but differ in detail. Furthermore, the 
shape of the massless apoT can be varied by adjusting 
S NP . At the same time, the massive resummed cross 
section is uniquely determined by our calculation. 

At sufficiently large Oh, both the massless and mas- 
sive resummed cross sections reduce to their correspond- 
ing finite-order counterparts. The massless cross sec- 
tion significantly overestimates the massive result near 
the threshold and at intermediate values of Q. For in- 
stance, at Q = 15 GeV (Fig. |(]a) the massless cross sec- 
tion is several times larger than the massive cross section 
in the whole range of Oh- In contrast, at Q = 50 GeV 
(Fig. |^b) the massless apo agrees well with the massive 
apo at Oh ^ 20° and overestimates the massive apo at 
Oh 20°. The massive apor is above the massless apo 
at 3° < H < 10° and below it at H < 3°. 

The presence of two critical angles (Oh ~ 3° and 
~ 10°) in apor can be qualitatively understood from 
the following argument. The rapid rise of apor over 
the massive ctfo begins when the small-^T logarithms 
ln m (qp/Q 2 ) become large — say, when qp is less than 
one tenth of Q 2 . Given that the Peterson fragmentation 
function peaks at about z ~ 0.95, and that Q = 50 GeV, 
M = 4.5 GeV, the condition ~ 0.1 Q 2 corresponds 
to Oh ~ 8°, which is close to the observed critical angle 
of 10°. Note that in that region q% > M 2 /z 2 . On the 
other hand, when q T is of order M 2 /z 2 , the growth of 
the logarithms ln(g^/Q 2 ) is inhibited by the non-zero 
mass term M 2 /z 2 in q%. The condition q T ~ M 2 /z 2 cor- 
responds to Oh ~ 2.5°, which is approximately where the 
mass-dependent cross section turns down. 
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VII. CONCLUSION AND OUTLOOK 

In this paper, we presented a method to describe po- 
lar angle distributions in heavy quark production in deep 
inelastic scattering. This method is realized in the sim- 
plified ACOT factorization scheme ||l7| , |25| and uses the 
impact parameter space (6-space) formalism |35|, ||(| 
to resum transverse momentum logarithms in the cur- 
rent fragmentation region. We discussed general fea- 
tures of this formalism and performed an explicit cal- 
culation of the resummed cross section for the O(a s ) 
flavor-excitation and 0{a 1 s ) flavor-creation subprocesses 
in bottom quark production. According to the numeri- 
cal results in Section VI, the multiple parton radiation 
effects in this process become important at Q > 15 GeV 
(or approximately at Q 2 /M 2 > 10). At Q = 50 GeV, 
the multiple parton radiation increases the inclusive cross 
section by about 25% as compared to the finite-order 
flavor-creation cross section. 

Many aspects of the resummation in the presence of the 
heavy quarks are similar to those in the massless resum- 
mation. In particular, it is possible to organize the calcu- 
lation in the massive case in a close analogy to the mass- 
less case by properly redefining the Lorentz invariants (in 
particular, by replacing the Lorentz-invariant transverse 
momentum qx in t he logarithm s by the rescaled trans- 
verse mass qr = V <$• + M 2 /z 2 ). The total resummed 
cross section is presented as a sum of the 6-space integral 
and the finite-order cross section <jfo, from which 
we subtract the asymptotic piece ctasy- Constructed 
in this way, the resummed cross section reduces to the 
finite-order cross section at Q w M and reproduces the 
massless resummed cross section at Q ^> M. 

At the same time, there are important differences be- 
tween the light- and heavy-hadron cases. For instance, 
the light hadron production is sensitive to the coher- 
ent QCD radiation with a wavelength of order 1/Aqcd, 
which is poorly known and has to be modeled by the 
phenomenological "nonperturbative Sudakov factor". In 
contrast, in the heavy-hadron case such long-distance ra- 
diation is suppressed by the large value of M. Hence, for 
a sufficiently heavy M, as in bottom quark production, 
the resummed cross section can be calculated without in- 
troducing the nonperturbative large-6 contributions. It 
will be interesting to test the hypothesis about the ab- 
sence of such long-distance contributions experimentally. 
Given the size of the differential cross sections obtained 



in Section |V| accurate tests of this approach will be fea- 
sible once the integrated luminosity of the HERA II run 
approaches 1 fb _1 . The same calculation can be done 
for charm production. However, in that case the region 
b > 1 GeV -1 is not as suppressed, and the nonperturba- 
tive Sudakov factor has to be included. 

Another important improvement in our calculation is 
more accurate treatment of threshold effects in fully dif- 
ferential cross sections. The accuracy in the threshold 
region is improved by introducing scaling variables ( |20| ) 
and ( ^l[ ) in finite-order and resummed differential cross 



sections. These scaling variables generalize the scaling 
variable proposed in Ref. for inclusive structure func- 
tions. They lead to stable theoretical predictions at all 
values of Q and agreement with the massless result at 
high energies. 

The extension of our calculation to higher orders is fea- 
sible in the near future, since many of its ingredients are 
already available in the literature ||| |7l). Further- 
more, in a forthcoming paper we will study the additional 
effects of threshold resummation [f| |7| ^ |7| ^ 
in DIS heavy-quark production, so that both transverse 
momentum and threshold logarithms are taken into ac- 
count. We conclude that the combined resummation of 
the mass-dependent logarithms \n(Al 2 /Q 2 ) and trans- 
verse momentum logarithms \n(q^/Q 2 ) is an important 
ingredient of the theoretical framework that aims at 
matching the growing precision of the world heavy-flavor 
data. 
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Appendix A: CALCULATION OF THE 
MASS-DEPENDENT C-FUNCTION 



In this Appendix, we derive the O(as) part of the 
function C l ^ G {x, fipb,bM). This is the only 0(as) term 

in the heavy-quark VF-term that explicitly depends on 
the heavy-quark mass M. This function appears in the 
factorized small-6 expression for the "6-dependent PDF" 

#'J A ( a;) fe,{m,},C 1 /C 2 ): 



hj A 



b,{ m q\^ ) = / ~T~~ 



^2/ J x t,a 



xC% a [x,n F b,bM;-^)f a/A [Z a , 



fJ-F 



(Al) 
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To perform this calculation, we consider a more ele- 
mentary form of Eq. (A.1), which represents the leading 
regions in Feynman graphs in the limit Q — ► oo. This 
elementary form can be found in Ref. J35|], where it was 
derived in the case of e + e~ hadroproduction. The func- 
tion ^"/a ( x , b, {m q }, C1/C2) is decomposed as 



— 1 a 

j/A 



<^2 



1 C 2 b 



xU(b) 1/2 ^J A [x,b,{m q },fi 



Ci 
C 2 



(A2) 



Here J#j denotes the "hard vertex", which contains highly 
off-shell subgraphs. U denotes soft subgraphs attached 
to Mj through gluon lines. ^p A (x, b, {m q }, C1/C2) 
consists of subgraphs corresponding to the propaga- 
tion of the incoming hadronic jet. The jet part 
^Jl A {x, b, {m q }, C1/C2) is related to the ^-dependent 
PDF ^p A (x,k T ,{m q },CA), defined as 



<?™ A (x,k T ,{m q } > tA) = Y l J2 J 



spin color * 



dy d 2 y T 
(2tt) 3 



ixp\y +ik T -y 1 



x(p A \ipj(0,y ,Yt)^4>-M\pa) 



(A3) 



in the frame where p A = {p A ,0,Or} , Pa = 
{xp A , M 2 /(2xp A ), k^} , and p\ — > 00. This definition 
is given in a gauge r\ ■ si = with rj 2 < 0. The kr- 
dependent PDF depends on the gauge through the pa- 
rameter^ = {PA ■ ri)/\r) 2 \. 

Let ^™i A (x,b,\mq\XA) be the 6-space transform of 
^Jj A (x, kr, {m,}, (a) taken in d dimensions: 



jd-2i ik T b 



^ A {x,b,U,{m q }) = J d d - 2 k T e" 
x&fi A (x,k T ,(A,{rn q }). (A4) 

Note that our definition t^Jj^x, b, ( A , {m q }) differs 
from the definition in Ref. Q by a factor (2ir) 2 ~ d . 
The jet part &%} A (x, b, {m q }, C 1 /C 2 ) is related to 
^jj A (x, b, {m q }, ( A ) in the limit ( A -> 00: 

&p A Lb, W,g) = Je 5 '^^ 



x&f} A (x,b,{m q },U)y 
where S'(b, (a', C1/C2) is a partial Sudakov factor, 

s'{b,U;C x /c 2 ) 



(A5) 



n An , 
c ^ dp, 



Ci/b f* 













-< — v 



PA* 



Figure 7: The Feynman diagram for the 0(cts/ir) contribu- 
tions f^j G {x,fi F /M) and lim^^^ 3?*™^ (x,b, M, ( A ). 



I ^4 — ) ^,^{a s {n)) 



/'• 



J? \ b\as 



Ci 
b 



C 2 



(A6) 
can be 



The definitions of the functions jjgr, Jf, and 
found in Ref. @. 

We now have all necessary ingredients for the calcula- 
tion of the 0(as/7r) function C l ^ G \x, bM). Setting 
j = h and A = G, and expanding Eqs. ( Al . A2 IAS!) , and 
3) in powers of as/V, we find 



C%fg(x,ii F b, bM) 



lim 



{&r$(x,b,M 



Xa)} 



(A7) 



where the superscript in parentheses denotes the order 
of as/ft. In the derivation of this equation, we used the 
following easily deducible equalities: 



,(0) 



1. 



s 



(o) 



4/i 0) (*) 



csi™(0) 
r(0) 



,m(0) 



(0) 
h/G 



f»; G (x) = s(x-i). 



(A8) 
0, (A9) 
(A10) 



The r.h.s. of Eq. ( |A7| ) can be calculated with 
the help of the definitions for fj^) G (x,fiF/M) and 

^/c\x, b, M, ( A ) in Eqs. (|) and respectively. 
A further simplification can be achieved by observing 
that at <D{a s /Tr) the limit rj 2 -> in &>H£(x, b, M, ( A ) 
can be safely taken before the limit (a — > 00, and, fur- 
thermore, for ?7 2 = the function (x, 6, M, £4) 
does not depend on (a- Correspondingly, both objects 
can be derived in the lightlike gauge from a single cut 
diagram shown in Fig. 0, where the double line corre- 
sponds to the factor r y + S(p A —p^ — k + )/2 in the case of 

fi% (x^f/M) and 7 +J(p+ - p+ - fc'+)e* k -' b /2 in the 
case of limbec &j$\x,b,M,U)- 

The difference between ]im^ A - too (x, 6, M, (a) = 

(x, b, M) and fh/Q (x, p,p/M) resides in the ex- 
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/ in(l) 

tra exponential factor e lk T b i n (x, b, M). Re- 

markably, this factor strongly affects the nature of 
^™, G (x,b,M). The loop integral over k T in 

fh/G ( x ' Mf/-&0 contains a UV singularity, which is regu- 
larized by an appropriate counterterm. In the ACOT 
scheme, the UV singularity is regularized in the MS 
scheme if /if > M, and by zero- momentum subtrac- 
tion if [If < M. The result for the heavy-quark PDF 
fff G (x,n F /M) is 



P, 



f/ G (x)ln^ F /M), ^ F >M; 



Ml 



0, fx F < M. 



(All) 



As expected, f^j G (x, fi F /M) exhibits the threshold be- 
havior at fx F = M. 

In contrast, the UV limit in the loop integral of 
^h/G ( x ' ^ is regularized by the oscillating expo- 
nent e ikT ' b . Since no UV singularity is present in 
^wc? (x, b, M), it does not depend on fi F and, there- 
fore, does not change at the threshold. It is given by 



0$\x,b,M) = Ptf G (x)K (bM) 

+ T R x(l — x)bMKi(bM). (Al2) 

Here KpjbM) and K\(bM) are the modified Bessel func- 
tions fl78], which satisfy the following useful properties: 



lim K (bM) 

bM^oa 



lim bMKi (bM) = 0; (A13) 

bM—*oo 



K (bM) -> - ]n(bM/b ) as bM -v 0; (A14) 
bM K\ (bM) -> 1 as bM -> 0. (A15) 



The "infrared-safe' 



*r /G ~'&b,M) 



part 
is obtained 
f§ G (x,n F /M) as in Eq. ©: 



x (l + Cl (bM)) 
+P { h %(x) (c (bM) 



by subtracting 



Trx(1 - x) 



(A16) 



-,m(l) 



= C I/G (x,^ F b,bM) 



/IF < A/ 

fj, F >M 



M ' 



(A17) 



In these equations, co(6M) and ci(bM) are the parts of 
if (&M) and bM Ki(bM) that vanish at 6M -> (cf. 
Eqs. (|Al3[fAT5|)): 



co(bM) = # (&M)+ln— ; (A18) 
ci(6M) = bMKi(bM) - 1. (A19) 



If /if is chosen to be of order 6q/6, no large logarithms 
appear in C h n , G (x, fx F b,bM) at b — * 0. At large Q, the 
small-6 region dominates the integration in Eq. (p"7|), so 
that C™, G (x, n F b,bM) effectively reduces to its massless 
expression |fl], |1: 



C$\x,r F b,bM) 
-P, 



6—0 



T R x(l - x) 



(A20) 



The above manipulations can be interpreted in the 
following way. At small b (b = bo/fi F < bo/M), we 

subtract from &™iq (x, b, M) its infrared-divergent part 



P^ G (i) ln(/ii?/M), which is then included and resummed 

in the heavy-quark PDF fh/G{ x i Mf /M). The convolu- 
tion of the resulting C-function with the PDF remains 
equal to ^T/q (%, b, M) up to higher-order corrections: 



h/G 



Enin 
L h/a 

a=h,G 



fa/c = & h /G (x, mf&, bM) + 0(a 2 s ). 

(A21) 



At large b (b > bo/M), the heavy-quark PDF fh/G is 
i dent ically equal to zero. To preserve the relationship 
(A21) below the threshold, one should include the above 
logarithmic term in the function C™,q (x, fipb, bM), as 
shown in Eq. (A17). The addition of an extra term 



P^/ G {x) In (p F /M) to C™£ J (x,fiFb,bM) at fip < M en- 
forces the smoothness of the form-factor W(b, Q, x, z) in 
the threshold region, which, in its turn, is needed to avoid 
unphysical oscillations of the cross section dajdq\. 



Appendix B: THE FINITE-ORDER CROSS 
SECTION 



This Appendix discusses the finite-order cross section 
da F o/d& that appears in the factorized hadronic cross 
section (ff7l). For the O(a s ) subprocess e + h — > e + h, 
this cross section is the same as in the massless case: 



da(e + h — ► e + h) 
d$ 



a Fi Ax($,ip) 



FO 



S, 



cA 



xe J 2 <5(q T )(5(l-x) ( 5(l-z), (Bl) 
where, in accordance with the notations of Ref. [B2|, 

Q 2 



Fi = 



AnS eA x 2 V 2 1 ' 
2 Q 2 ' 



(B2) 



The contribution of the gluon-photon fusion channel is 
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dd(e + G -> e + h + h) 



FO 



a F t a s 2 
1 — o eg' 5 



1 

- - 1 

X 



1 

- - 1 



eft \ x{\ — x) 



:T R %{% Q 2 ,z, g 2 , M 2 )A p (ijj, tp), 
P =i 



(B3) 



where A p (ip, ip) denote orthonormal functions of the lep- 
tonic azimuthal angle tp and boost parameter ip given in 
Eq. (||). The structure functions V p (x,Q 2 ,z,q^,M 2 ) are 
calculated to be 



Vi. 



1 



x 2 cfe 



1 - 2xz + 2x z 



2c2 



M 2 x 2 



2z 1 
x Q 



— -p^ (§XZ — X — z) 



[ 4^z 2 + 2-8z + 8z 2 -4- 



M 



, (B4) 



V 3 



M 2 1 

4— — +4ki (-l + z)£, (B5) 
W" It 



■1 + 2 1 



9t 
Q 2 



xz 



q T 



4«iz(-l + 2S)^-, 



In Eqs. (B4-B7), 



M 2 (l -a) 



(B6) 
(B7) 

(B8) 



Appendix C: KINEMATIC AL CORRECTION 

In this Appendix, we derive the kinematical corrections 
(p0|) and (pll) that are introduced in the flavor-excitation 



contributions to <tfo, as well as in and oasy- Let 
us first consider the 0(as) cross section for the photon- 
gluon fusion, which we write as 



da{e + A^e + H + X) 



6 J £a 



Dh/Ii (6>) fa/ A (£a) 



*-i]-^l«»> 



(Cl) 



Here /?($) includes all terms in the parton-level cross sec- 



tion (d3 / d<&)FO except for the 5— function (cf. Eq. (B3)): 



T R ]T V p {x, Q 2 , z, q 2 Tl M 2 )A p {ip, <p). (C2) 
P =i 



The (S-function in Eq. (Cl) can be reorganized as 



'H(H-§' 



Z Q 2 



t¥ 4 - 4Af 2 (g 2 + W 2 



where 



W 2 ± ^W 74 - 4M 2 (g 2 + W 72 ) 
1 2M 2 : 



(C3) 



(C4) 



and W 2 = Q 2 (1 — x) /a;. We see that the mass- 
dependent phase space element contains two (^-functions 
5(£b — ) an d — £(T)i which can be used to integrate 
out the dependence on ^ in Eq. (Cl). 

It can be further shown that in the massless limit the 
solutions £f, = and £b = correspond to the heavy 
quarks produced in the current and target fragmentation 
regions, respectively. When M — > 0, the relationship 
simplifies to 



[ i-i 

X 



zQ 2 
W 2 



z J Q 2 ) 



where 



?6 




0(M 2 ) 



(C5) 



(C6) 



(C7) 



In this limit, the solution diverges (and, therefore, 
will not contribute) unless z is identically zero. How- 
ever, according to Eq. (|39|) and the last paragraph in 
Subsection B of Section |V|, at z = the observed 
final-state hadron appears among remnants of the tar- 
get {Oh ~ 180° in the 7* A cm. frame), i.e., away from 
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the region of our primary interest (small and intermedi- 
ate Oh)- Hence, in the limit Oh — > all dominant loga- 
rithmic contributions as well as their all-order sums (the 
flavor-excitation cross section and W-tevva) arise only 
from terms proportional to 6(£b — £&~ ) • The contributions 

I 



proportional to S(^b - 
region are suppressed. 



£k ) in the current fragmentation 



The integration over £b with the help of Eq. ( |C3|) l eads 
to the following expression for the cross section (CI): 



(do{e + e + H + X)\ f C - 
\ dxdQ 2 dzdq^dip ) , G J^ m , 



Q 2 



W 4 -AM 2 [q 2 +W 



I 



(C8) 



Here the lower and upper integration limits £™" and 
£,™ ax are determined by demanding the argument of the 
square root in Eq. (C£) be non-negative and £& < 1; that 
is, 




(C9) 



(CIO) 

for £b = £b ■ We see that, according to the exact kinemat- 
ics of heavy flavor production, the heavy quark pairs are 
produced only when the light-cone momentum fraction 
£ a is not less than £™ n (where £™ m > x) and not more 
than C iax (where ^ ax < 1). The exact values of £™ n 
and £,a iax are different for the branches with = ^ and 

Turning now to the flavor-excitation contributions 7* + 
h — > h + X, we find that in those the integration over £ a 
a priori covers the whole range x < £ a < I. Indeed, in 
those contributions the heavy antiquark in the remnants 
of the incident hadron is ignored, so that the reaction 
can go at a lower cm. energy W than it is allowed by the 
exact kinematics. Since the PDF's grow rapidly at small 
x, the naively calculated total cross section otot tends to 
contain large contributions from the unphysical region of 
small x and disa gree with the data. To fix this problem, 
we use Eq. ( CIO ) to derive the following scaling variable 
in the finite-order flavor-excitation contributions: 



Xh = x 1 + 



1 



Ar 



z(l-z) 



(Cll) 



This variable takes into account the fact that the incom- 
ing heavy quark in the flavor-excitation process appears 



from the contributions with = in the flavor-creation 
process, and that the transverse momentum zqr of this 
quark in the finite-order cross section is identically zero. 

Similarly, we notice that the W-term and its finite- 
order expansion oasy contain the "6-dependent PDF's" 

^KIA ( X A M, Ci/C 2 ), which correspond to the incom- 
ing heavy quarks wi th a non-zero transverse momentum. 
According to Eq. ( C1C| ), the available phase space in 
the longitudinal direction is a decreasing function of the 
transverse momentum zqt, and it is desirable to imple- 
ment this phase-space reduction to improve the cancella- 
tion between o^ and oasy at large qr- In our calcula- 
tion, this feature is implemented by evaluating o^ and 
oasy at the scaling variable 



Xh 



= x\l 



M 2 



z(l 



Q 



(C12) 



which immediately follows from Eq. ( CIO ). 

Des pite the a ppare nt complexity of the scaling vari- 
ables ( Cll ) and ( C12 ) , they satisfy the following impor- 
tant properties: 

1. They are straightforwardly derived from the ex- 
act kinematical constraints on the variable £ a in 
Eqs. @ and jC10|) . 



They remove contributions from unphysical values 
of x at all values of Q and qr, thus leading to better 
agreement with the data. 



3. In the limit Q 2 > M 2 



th e variable Xh m a FO re- 
duces to x (cf. Eq. ( Cll )), so that the standard 
factorization for the massless finite-order cross sec- 
tions is reproduced. 



4. In the limit Q 2 > M : 



and Q 2 3> q 2 -,, the variable 



Xh m a w and asy reduces to x (cf. Eq. (|C12|)), 



so that the exact resummed cross section is repro- 
duced. 
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Finally, consider the integration of the cross section ( |Cq ) 
over z, q\, and ip to obtain the 0(as) 7*G contribution 
to an inclusive DIS function F(x,Q 2 ). We find that 



F(x,Q 2 



1 d£a r (i) ( X \l F M 

x fG < A (*»{%})' (C13) 



where the lower limit of the integral over £ a is given by 

4M 2 ' 



£a = X 1 



Q 2 



(C14) 



£ Q can be easily found from Eqs. (C9) and ( |C10| ), given 
that q T > 0, < z < 1, and z(l - z) < 1/4 in the 
interval < z < 1. Since in the 7*G contribution the 
integration over £ a is constrained from below by £ a > x, 
it makes sense to implement a similar constraint in the 
flavor-excitation contributions by introducing the scal- 
ing variable Xh = x(l + AM 2 /Q 2 ). This variable is pre- 
cisely the one that appears in the recent version of the 
ACOT scheme with the optimized treatment of the in- 
clusive structure functions in the threshold region [p9[ . 
Our scaling variables extend the idea of Ref. j29| to the 
semi-inclusive and resummed cross sections. 



for both solutions = and £& = £, b ■ This value of 
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